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We study a diffusive superconductor-ferromagnet-superconductor (SFS) junction with in-plane 
ferromagnetic domains. Close to the superconducting transition temperature, we describe the prax- 
es ■ imity effect in the junction with the linearized Usadel equations. We find that properties of such a 
junction depend on the size of the domains relative to the magnetic coherence length. In the case 
of large domains, the junction exhibits transitions to the ty state, similarly to a single-domain SFS 
junction. In the case of small domains, the magnetization effectively averages out, and the junc- 
tion is always in the zero state, similarly to a superconductor-normal metal-superconductor (SNS) 
junction. In both those regimes, the influence of domain walls may be approximately described as 
h— 5 ■ an effective spin-flip scattering. We also study the inhomogeneous distribution of the local current 
' density in the junction. Close to the 0-7r transitions, the directions of the critical current may be 
I i opposite in the vicinity of the domain wall and in the middle of the domains. 

£h \ I. INTRODUCTION 

4h ■ 

Substantial progress has been made recently in understanding the physical properties of nanoscopic layered struc- 
tures composed of superconducting (S) and ferromagnetic (F) materials (for a review, see Refs. 1,2). In such systems, 
the coexistence of superconducting and magnetic correlations may lead to a variety of interesting physical effects. The 
exchange splitting of the Fermi level 3 in the ferromagnet breaks the spin degeneracy, and Andreev-reflected Cooper 
pairs acquire a finite momentum, which produces oscillations of the Cooper-pair wave function. By tuning the geo- 
metric and electronic parameters, one can realize SFS junctions in which the superconducting order parameter is of 
opposite sign in the two S electrodes. 4 Such a tt state manifests itself in a reversal of the sign of the critical current 
. as the thickness of the ferromagnets is varied. 1 The characteristic thickness at which the first 0-7r transition occurs is 
of the order of the magnetic coherence length In the diffusive limit, which is realized in most experimentally fab- 
ricated SF heterostructures, £/, is given by y/ D/h where D denotes the diffusion constant and h is the ferromagnetic 
exchange energy. Therefore, the experimental observation of such 0-7r transitions in nanoscale devices requires a low 
exchange energy h. This stringent condition was achieved using weak ferromagnetic CuNi or PdNi alloys. 5-8 

The physics of single-domain SFS junctions (including the effect of spin-flip 9-11 and spin-orbit scattering 3,12 ' 13 ) 
is now well understood. However, in some experiments the 0-7r transition points may deviate from standard 
predictions 8,14 or even be absent. 15,16 There is no consensus on the interpretation of such deviations. It may be 
attributed either to the presence of a magnetically dead layer at the interface between the superconductor and the 
ferromagnet, 8,17 or to a domain structure or inhomogeneities in the ferromagnetic layer. The domain structure cru- 
cially depends on the nature of the ferromagnet: strong ferromagnets consist of well-defined magnetic domains whose 
spatial extension may be reduced by the proximity effect 18-21 . In weakly ferromagnetic alloys, on the other hand, the 
magnetization may fluctuate on short length scales without forming domains. 22 

Theoretically, SFS junctions with inhomogeneous magnetization have been studied recently in different setups. 23-28 
• • . However, in most works (except for Refs. 27,28) only domains along the junction were studied (quasi-one-dimensional 
geometry), while in the experimental realizations of SFS junctions with thin F layers the domain structure is more 
likely to form in the plane of the F film. 

Motivated by the experimental progress on ir junctions, we study a model of a diffusive SFS junction with in-plane 
domains (so that the domain walls are orthogonal to the S and F layers, see Fig. 1). This geometry has been studied 
previously by Volkov and Anishchanka within the macroscopic approach of London equations. 28 Our model is different 
from the one studied in Ref. 27: in that work, the Neel domain walls are considered, and the junction is brought to the 
regime with only the long-range triplet component contributing to the Josephson current. In our model, the domain 
walls are taken to be sharp, and no long-range triplet component appears for domains with antiparallel magnetization. 

The domain structure introduces an additional length scale: the domain size a. As one can expect, we find that the 
effect of inhomogeneous magnetization depends strongly on the relative magnitude of a and In the limit of small 
domains, a <C the exchange field effectively averages out, and the critical current of a single nonmagnetic SNS 
junction is retrieved. In the opposite limit of large domains a 3> £/u the influence of domain walls is localized to their 
vicinity and produces only a small correction to the current of a single-domain SFS junction. Between those limits, 
the supercurrent shows either a damped oscillatory behavior as a function of the junction thickness (for large domains 
a > a c w 0.83 £,h), or a monotonic exponential decay (for smaller domains a < a c ). In the former case, the multidomain 
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FIG. 1: SFS junction with in-plane magnetic domains 



junction may be compared to a single-domain SFS trilayer with spin-flip scattering 12 and a renormalized exchange 
field, whereas in the latter case the junction behaves like a SNS junction with spin-flip scattering. The effective 
parameters are determined analytically in both limits of small and large domains. We also study the inhomogeneous 
distribution of the current density and conjecture that at low temperatures such SFS junctions with domains may 
realize the intermediate <p phase proposed by Buzdin. 29 

The paper is organized as follows. In Section II we compute the superconducting Green functions and the Josephson- 
current density for the multidomain SFS junction. Section III is devoted to the analysis of the total Joscphson current. 
In Section IV we discuss the spatial distribution of the current density. Finally, in Section V we summarize our 
conclusions. 



II. A MODEL FOR THE MULTIDOMAIN SFS JUNCTION 



We assume that the ferromagnetic layer is strongly disordered, and the motion of electrons is diffusive. In this regime, 
the Green functions are given by the solutions to the Usadel equations. 30 To simplify the calculations, we further 
assume that the junction is close to the superconducting critical temperature T c . In this case, the superconducting 
correlations are weak so that the Usadel equations can be linearized, and the current-phase relation is sinusoidal 31 

J = J c smip, (1) 

where ip = 2% is the superconducting phase difference across the junction and J c is the critical current. The sign of 
J c determines if the junction is in the "zero" phase or in the 'V phase. 

In this paper, we consider a SFS junction with in-plane ferromagnetic domains of opposite magnetization. We 
introduce a coordinate system with the F layer in the yz plane (Fig. 1). The x axis is directed along the junction, 
and the SF interfaces correspond to the coordinates x = 0, d. The domain walls are taken to be normal to the y axis. 
The origin of the y axis is chosen at the interface between two domains. The system is invariant under translation 
along the z axis. Our further calculations will be equally applicable to either the system with two domains of width 
a (see Fig. 1) or the 2a-periodic multidomain case (the same setup periodically repeated in the y direction). 

The (nonlinear) Usadel equation in the ferromagnetic layer takes the form (we follow the conventions used in Ref. 
32) 

(gWg) - lo [f 3 «J , g] - i [f 3 (h • 5) , g) = 0. (2) 

where D denotes the diffusion constant and the system of units with ft=fcs=/iB = lis chosen. The Green function 
g is a matrix in the Nambu ® spin space, f a and a a denote the Pauli matrices respectively in Nambu (particle-hole) 
and spin space, to = (2n + 1) ttT are the Matsubara frequencies and h is the exchange field in the ferromagnet. The 
Usadel equation is supplemented with the normalization condition for the quasiclassical Green function 



g 2 = 1 = t <7 . 



(3) 
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For simplicity, we assume that the superconductors are much less disordered than the ferromagnet, and then we can 
impose the rigid boundary conditions at the SF interfaces, 

1 ( u Ae ±l *\ . ... 

9= 1 ,, , A9 1 ) ® CT o, (4) 

Nambu 



Vlu 2 + A 2 V -Ae™ -w 



where A denotes the superconducting order parameter, and the different signs refer respectively to the boundary 
conditions at x = and x = d. At the (transparent) interface between the two ferromagnetic domains, we impose the 
continuity of the Green functions and of their derivatives. 

Close to the critical temperature T c , we linearize the Usadel equations (2), (3) around the solution for the normal 
metal state g = T3<JoSgn(a;). The linearized Green function then takes the form 

/ er sgn(w) f a a a \ . . 

9 \ -fi* a -<7 sgn(o;) ) ' (5) 

where the scalar f (respectively /q) and vector f (respectively ft) components of the anomalous Green functions 
obey the linear equations 

92 d2 ^ -i(t) \2 ,i(t) 



dx 2 dy 2 



f^>-Xif^>=0, (6) 



with 

A± 



\u\ T ift-sgn(w) 1 ] 



D 



(7) 



The projections of the anomalous Green function along the direction of the exchange field ("parallel" components) 

are defined as f±^ (x,y) = ± f ^) • e z (we assume that the ferromagnetic exchange field h is aligned in the direction 
e z , see Fig. 1). Note that there is no perpendicular ("long-range triplet" 33 ) component of the vector part of the Green 
function, since the magnetizations of the domains are collinear. We have used the invariance under translation along 
the z direction. The superscript 1 refers to domains with field along e z and in the following we will use the superscript 

2(f) 

2 for domains with the field along — e z . Similar equations hold for f± with X± <-> X^. 
It is convenient to write the solutions to those equations in the form 

f2 n '\x,y) = f^tlW+S^^y) (8) 

where / ± ' Bulk are the solutions of Eq. (6) for a single-domain SFS junction with the magnetization along e z (respec- 
tively — e z ). Since the equations (6) are linear, the correction S(x, y) is also a solution to the same equations with the 
boundary conditions 

d y 5]l 2 (x,y = T a) = 0, (9) 

S]: 2 (x = {0,d},y) = 0, (10) 

Si(x, y = 0)-5 2 ± (x,y = 0) = Af ±Bulk (x), (11) 

[dy] 5i (x, y = 0) - [8y] 5 2 ± (x, y = 0) = 0. (12) 

Here A/±Buik = .f±Buik( x ) — /±Buik( x ) i s tn c difference of the bulk Green functions in the two domains. For the two- 
domain junction, the first condition imposes zero current at the interface with vacuum, the second condition ensures 
the continuity of the Green functions at the SF interfaces. Finally, the last two conditions reflect the continuity of the 
Green function and its derivatives at the interface between the two domains. It can be easily shown from symmetry 
considerations that this set of boundary conditions can also be applied to a periodic multidomain SFS junction with 
domains of width 2a. 

The condition (10) allows us to express d± 2 in the form of the Fourier series 

4' 2 = Esin(^)^( y )- (13) 



For each n we solve 



dX±=7^ ± (14) 
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with 

7™± = \/(^) 2 + A 2 ± . (15) 



To obtain the equation for A^± one needs to substitute ■y n ± <-> 7 n =p- We can solve those equations for each Fourier 
component n with the boundary conditions provided by (9), (11) and (12). The solution is given by 

ii = Af sin (HUE.) g^cosh 7 „ ± ( y + q) 7KT tanh 7KT a / 1 _ 1 \ 

l w l tZi cosh7„±a 7„ T tanh7„ T a + 7„± tanh7„±a \7^- 7± / 

In the second domain, the correction 5± is given by the same formula with the replacement of y, j± by —y, j T . The 
bulk Green functions are given by 26 



f 1 - A 
/±Buik - rr 



sinhA±a; _ ix sinhA±(d — x) f 
sinhA±d sinhA±<i 



(17) 



and j±Buik = /ipBuik- Finally, note that / Bu i k an d ^ are given by the same expressions (16), (17) with the replacement 
of X by -x- 

The last step will be to compute the Josephson current density using the formula 2 

oo ^ 

J = ieN(0)DnT E -Tr (T 3 a gVg) , (18) 

UJ — — oo 

where iV(0) is the density of states in the normal metal phase (per one spin projection) and the trace has to be taken 
over the Nambu and spin indices. The current density can be explicitly rewritten for the linearized g 



J = -ieN(0)DnT ]T 



- ± 2 



(19) 



The symmetry of translation along the z direction implies that the current remains in the xy plane. Using the 
expression for the Green functions (16) and (17), we can obtain a general expression for the current density (which is too 
cumbersome to be reproduced here). This expression involves two contributions. The first one is produced exclusively 
by the bulk Green functions (17) and corresponds to a homogeneous ferromagnetic interlayer. The second contribution 
is due to the correction (16) and reflects the influence of the domain structure. The current resulting from this 
contribution is not uniform in space. The characteristic decay scale of this correction as a function of the distance from 

the domain interface is given by 3? (^f^j ~ mm (£T, 6i, eQ, where £t = \/D/2itT c and £/, = ^D/h are the thermal 

and magnetic coherence lengths, respectively. Far from the interface between the domains (y 3> min(£-r, £h, d)), the 
correction (16) vanishes and we recover locally the single-domain SFS current. Thus we expect the properties of the 
junction to be very different in the two limits of small [a <C min(£r, 6u d)] and large [a >• min(£r, d)] domains. 



III. CRITICAL CURRENT 



Experimentally, in SFS hybrid junctions, the measurable quantity is the total current flowing through the junction, 
that is along the x-axis. Since V • J — 0, the total current is conserved along the x direction. We can therefore 
compute it at x = 0, and we find 



Jc 

h 



A 2 



.^>o 



uj 2 sinhA+d 
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167^ 

ad 2 C, 



A 2 

U1>0 



E 

71=1 
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(7«+7n-) 3 In- coth7„+a + 7„+ coth7„_a 



(20) 



with 



4ejV(0)£>5'7rr 
d ' 



(21) 



and S the area of the junction. The first term is the critical current for a single-domain SFS junction with a damped 
oscillatory dependence on the F-layer thickness (for a review, see Rcfs. 1, 34 and 35). It can be either positive 
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(zero state of the junction) or negative (ir state). The second term reflects the influence of the domain structure. 
The critical current (20) depends on the three dimensionless parameters: a/^ h , d/£ h , and £t/6i- For some values of 
the parameters, the critical current (20) computed numerically is plotted in Fig. 2. Depending on the values of the 
parameters, it shows either an exponential decay or an exponential decay with oscillations, as a function of d. 

Note that in most experimental situations £t 3> £h, because the ferromagnetic exchange energy exceeds by far the 
superconducting critical temperature. In the following we will refer to this situation as the high-field limit. In this 
limit, the summation over u in Eq. (20) can be performed analytically E w>0 ^ 2 /uj 2 — A 2 /(8T 2 )], and the deviation 
SJ C from the critical current of a single-domain SFS junction [the second term in Eq. (20)] is expressed in terms of 
the reduced variables n* = ^2 an d a * = 

la "2T 2 a* £j („4 + n * 4 ) 3 / 2 9? [Vn 2 + in* 2 coth (a*Vn 2 - in* 2 )} ' 

In the limit of large d, the asymptotic behavior of this expression may be estimated as an integral (in the variable 

z = n/n*) 



oo 

S_Jc = _^_ f dz e*™ z z 2 

h " 2T 2 a* J Z (z 4 + l) 3 / 2 [Vi 2_ Hcoth(a*n*v / i 2 ^) + Vz^^i coth (a* n* Vz^Ti)] ' 



which is, in turn, determined to the exponential precision by the singularities of the integrand in the complex plane. 
Remarkably, the contribution from the poles at (iz) 1 / 2 cancels exactly the first term (single-domain contribution) of 
Eq. (20). For sufficiently large d, to the exponential precision, the critical current is then given by 

Jc«e-", A = -^, (24) 

where zq is the singularity of the integrand with the smallest positive imaginary part. Note that zq is now a function 
of one dimensionless parameter a = a*n* — ^/2a/^h- 

By analogy with a single-domain SFS junction with spin-flip scattering, the real and imaginary parts of A 2 may be 
interpreted as an effective magnetic field and an effective spin-flip rate 38 , 



(off) 2 T D ' Sfc V /i( eff ) ' 



Therefore the effective field and spin-flip rate can be found as 



■ sf - 2 

In the following, we discuss the limits of large and small domain sizes. 



= *>(*g) , rif = _/W). (26) 



12 3a/i 



FIG. 3: Phase diagram of the junction in the high-exchange-field limit. Here a represents the width of the domains and d is 
the length of the junction. 



A. Limit of large domains : a>^ 

We consider the limit of large domains, a>((„ with the assumption of the strong exchange field, £t 3> In this 
regime, the damped oscillations of the critical current at large d are determined by the solutions to the equation 



V z 2 + i coth(av z 2 — i) + V ' z 2 — i coth(a-\/ z 2 + i) = 



(27) 



with the smallest positive imaginary part. At a = ^/2a/^h 3> 1, one of the arguments of coth(aV z 2 ± 1) must be 

close to ±i7r/2. Expanding around this point, we obtain Zq — i — + ^ 4 ^ i 7r — \- . . . This translates into the reduced 
effective field 



h (eS) 



8V2 



(28) 



and the effective spin-flip rate 



,(eff) 



7T J 

Tg 



tt 2 D 
16a 2 



(29) 



Thus, to the leading order in (^/t/a), the effect of domain walls reduces to an effective spin- flip rate, which increases 
the period of 0-7r transitions as a function of d and simultaneously decreases the overall decay length of the critical 
current (see Fig. 2, dashed line, for an illustration). 



B. Limit of small domains : a <C £,h,d, £t 



In the limit of small domains a <C £t, we can calculate a perturbative correction to the critical current by 
expanding (20) in a. To the lowest order in a, we obtain (without assuming the high-field limit), 



Jc 



JSNS 



2a 2 d 2 v A 2 



Xod cosh Xod — sinh Xod 



Aodsinh Xod 



(30) 



where Aq 



_ 2M 



-jy does not contain the exchange energy 



h, and Jsns = Jc(h = 0). This expression reveals 
the multidomain SFS junction behaves like a SNS junction: the exchange field is averaged out 
Note also that the correction arising from a finite domain width is always negative: 



that in the limit a 
when the domain width is small 
the amplitude of the current is decreased compared to the SNS case 
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A more accurate approximation may be obtained in the high-field limit £t S> £,h by the asymptotic estimate of the 
oscillating sum described earlier in this Section. To the second order in a, the solution to the equation (27) is given 

2 

by zq = — which translates into 



, , \ 2 

1 I a 



This expression for T^f' agrees with the general estimate for the effective spin-flip rate obtained by Ivanov and 

Fominov 32 for SF structures with inhomogeneous magnetization. 

Note that for sufficiently small a, the equation (27) has a solution with real Zq corresponding to a pure decay 
(without oscillations) of the critical current. The dependence of the critical current on d is then purely decaying, 
without 0-7r oscillations (Fig. 2, dotted line). 



C. 0— 7r phase diagram 



Between the two regimes of small and large domains, there is a phase transition as a function of a/£h corresponding 
to a bifurcation of the real solution Zq to Eq. (27) at smaller a to complex solutions at larger a. For a/£h smaller than 
the critical value, the critical current decays as a function of d without oscillations (always in the phase). For a/£h 
larger than the critical value, the dependence on d is damped oscillatory, qualitatively similar to a single-domain SFS 
junction. 

Numerically, we find the critical value a c /£h « 0.83. The full 0-7r phase diagram in the high- field limit is plotted in 
Fig. 3. Periodic 0-7T transitions (as a function of d) above a c /^h and zero phase below a c /^h illustrate our discussion. 
The absence of the 0-7r transitions in the case of small domains may explain why in some experimental SFS junctions 
the 7r state is absent. 15,16 

For completeness, in Fig. 4 we also plot the locus of solutions Zq to Eq. (27) in the complex plane for all values of 
a (in the units T^F /h and ft,( cff ' /h). The corresponding real and imaginary parts of A determining the d dependence 
of the critical current (24) are plotted in the inset. 



IV. LOCAL CURRENT DENSITY 



Since the system does not have a translational symmetry along the y direction, the Josephson current forms a 
nontrivial pattern in the x-y plane. In Fig. 5 we present plots of the current density (proportional to sin ip) at two 
different points of the phase diagram: in the zero phase and in the 7r phase. 

Those inhomogeneous patterns may be qualitatively understood on the basis of interpreting the domain walls as 
producing an effective spin-flip scattering. Different regions of the ferromagnet may be attributed different effective 
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FIG. 5: Josephson-current density in the two-domain SFS junction in the zero and n phases. The domain size is a = 1.6 £h <C £t, 
and d is taken to be below (left) and above (right) the first 0-7r transition. 




FIG. 6: Current lines close to the 0-7r transition. The domain size is a — 1.6 £h <C £t, and d is taken to be close to the first 
0-7r transition. 



spin-flip rates, depending on their distance from the domain wall. The effective spin-flip processes renormalize the 
decay coefficient A in (24) and, therefore, different parts of the junction experience 0-7T transitions at different values 
of d. This can be clearly seen in Fig. 6 depicting the current density near a 0-7T transition. While the neighborhood 
of the domain wall is in the phase, the region near the free boundaries (at y = ±a) are in the it phase. This 
situation resembles a model studied by Buzdin et al. 29 : a system of alternating zero and ir junctions. In that work, an 
intermediate equilibrium phase difference was predicted, depending on the ratio between the junction widths and the 
magnetic coherence length. Even though our model cannot lead to such a ^-junction (we consider linearized Usadel 
equations and therefore obtain a purely sinusoidal current-phase relation with only two possible equilibrium phases 
or 7r), at low temperatures such a SFS system with domains could possibly produce a estate. 



V. SUMMARY 



In this work we consider a Joscphson SFS junction consisting of domains with opposite magnetization connected 
"in parallel". As a function of the junction thickness, the critical current may exhibit either a decaying oscillating 
or a purely decaying behavior, depending on the domain width. The effect of domain walls in this geometry may 
be approximated as an effective spin-flip scattering, together with a renormalization of the effective magnetic field. 
This behavior is different from that in SFF'S junctions with the domains connected "in series" studied in Ref. 26. In 
that SFF'S setup, the domain structure lead to a gradual reduction of the n phase (at a non-parallel configuration 
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of the two domains), so that the relative fraction of the zero phase increases as a function of the mismatch in the 
magnetization directions. In the present work, however, we do not consider the case of an arbitrary angle between 
the two magnetizations, because of the complexity of the problem. 

We expect that in a realistic geometry of domains both effects of the spin-flip scattering and of the reduction of 
the 7r phase take place simultaneously, and our findings from this work and from Ref. 26 may help to qualitatively 
describe the O-tt phase diagram of real SFS junctions with inhomogeneous ferromagnets. 
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